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The photon emissivity from the bremsstrahlung process e~e~ — > e~e~7 occuring in the electrosphere at 
the bare surface of a strange quark star is calculated. For surface temperatures T < lO^K, the photon flux 
exceeds that of e^e~ pairs that are produced via the Schwinger mechanism in the presence of a strong electric 
field that binds electrons to the surface of the quark star. The average energy of photons emitted from the 
bremsstrahlung process can be 0.5 MeV or more, which is larger than that in e'^e' pair annihilation. The 
observation of this distinctive photon spectrum would constitute an unmistakable signature of a strange quark 
star and shed light on color superconductivity at stellar densities. 

PACS: 26.60.+C, 95.30.Cq, 97.60.Jd 



I. INTRODUCTION 



Recently, it has been pointed out that thermal emission from the bare surface of a strange quark star, due to both photons and 
pair production, can produce luminosities well above the Eddington limit (~ 10^^ erg sec~^) for extended periods of time, 
fi^m about a day to decades, depending on the superconducting phase of quark matter [1]. The spectrum of emitted photons 
Is^ignificantly different from that of a normal cooling neutron star (30 < E/keV < 500 instead of 0.1 < E/keY < 2.5). This 
SMinctive spectrum and temperature evolution, if observed, would constitute an unmistakable detection of a strange quark star 
and shed light on color superconductivity at stellar densities. The predicted characteristics are well within the capabilities of the 
&TEGRAL satellite [2] launched toward the end of 2002. 
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strange quark matter (SQM) star differs from a normal neutron star (NS) in several important ways, which include 



3 A SQM star is self-bound [3-7] while a NS requires gravitational forces for its binding. As a consequence, the pressure vanishes 
■yanishing baryon density in a NS with a surface of normal matter. The interior of the NS, however, may contain any or a 
cotnbination of non-nucleonic matter such as hyperons, pion or kaon condensates, or quark matter [8,9]. In contrast, a self-bound 
star is made entirely of quark matter up to its bare surface and the pressure vanishes at a finite but supra-nuclear baryon 
diMsity. In the context of the MIT bag model with first order corrections due to gluon exchange, the baryon density at which 
t^ pressure vanishes is [10] 

^b{p^Q)=[^A 1--^ , (1) 



\37r2/-^/ \ TT 



. wfiere B is the bag constant and ttc = 5c /(^tt) is the quark-gluon coupling constant. For typical values of B and ac [5], the baryon 
density at vanishing pressure is about 2 to 3no, where no = 0.16 fm^^ is the saturation density of nuclear matter. The density in 
Eq. (1) is not significantly affected by the finite strange quark mass [10] or by pairing gaps in the quark phase [11]. 

(2) The mass versus radius relations of NS and SQM stars differ significantly, although in the range of masses (1 < M/Mq < 2) 
observed to date [12], the calculated radii are similar {R ~ 10 to 15 km). This makes it difficult to distinguish these two classes 
on the basis of their gross physical properties alone. 

The light curves of NS and SQM stars, determined by surface photon emission, can however be very different under certain 
conditions. Central to the discussion of photon emission from a SQM star is the physical size and composition of its surface layer. 
If quark matter is in the Color-Flavor-Locked (CFL) superconducting phase (in which all three flavors take part in pairing) up to 
the surface, no electrons are required or admitted [13-15]. Surface photon emission from the CFL phase has been found to saturate 
the blackbody limit at early times [16]. In contrast, charge neutrality cannot be satisfied by three flavors of quarks alone in the 
two-flavor superconducting phase (2SC). In this case, an electron concentration Ue/nB of about 10^"* to 10^'^ is required to achieve 
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charge neutrality. Normally, the bare surface of a quark star is an inefficient emitter of thermal radiation below temperatures 
T ~ lO^^K because the plasma frequency of quark matter is high, of order 20 MeV [17]. However, in the 2SC phase, strong 
electric fields can arise because electrons in the surface layer are bound by electrostatic interaction to quark matter. This layer of 
electrons, which we will call the "electrosphere" , is typically 10^ fm thick [17] and can have radial electric fields whose magnitude 
(~ 5 X 10^^ V cm~^) exceeds the critical value required for electron-positron pair production from the QED vacuum. The critical 
value for the electric field was estimated by Schwinger to be [18] 

2 3 

E^, = ~ 1.3 X IQi'^ V cm-i , (2) 

en 

where me and e are the electron's mass and charge, respectively. For a homogeneous electric field E ^ E„, the vacuum pair- 
creation rate per unit volume is given by [18] 

Usov [19] first exploited the possibility of thermal emission of hard photons from the annihilation of e+e" pairs created by super- 
critical electric fields in degenerate matter and demonstrated that the associated photon emissivity is large enough to be observable. 
Subsequently, Page and Usov [1] have shown that the bare surface of a quark star under super-critical electric fields leads to a 
light curve that is significantly different from those produced by either quark matter surrounded by a layer of normal matter, or 
by surfaces of compact objects made entirely of normal matter. Their principal conclusion was that hard photon emission from 
e+e" annihilation can dominate significantly over the blackbody spectrum at lower luminosities, and that the variation in the 
photon spectrum with luminosity could be used as an observational signature of a young bare strange quark star. Furthermore, 
calculations of light curves from quark matter in its superconducting phases (CFL or 2SC) have shown that it may be possible to 
set bounds on quark pairing gaps provided they are in the range 0.5 MeV to a few MeV [1]. 

Our objectives in this paper are to 

(1) point out that, depending on the temperature, the bremsstrahlung process e~e~ e~e~7 in the electrosphere can lead to 
photon luminosities that are significantly larger than the e+e~ pair production luminosities in the presence of strong electric fields, 

(2) present a calculation of the corresponding emissivity including minimal screening effects, and investigate the role of the Landau- 
Pomeranchuk-Migdal (LPM) effect, and 

(3) compare the photon emissivity from the er er e~e~7 bremsstrahlung process to other processes occurring in the thin 
electrosphere at the surface of the star. Examples of such processes include e+e~ pair production which occurs due to the 
presence of strong electric fields, and equilibrium and non-thermal bremsstrahlung radiation from quark-quark collisions in the 
uppermost layer of quark matter. 

We wish to emphasize that our discussion applies only to the case in which a charge neutralizing surface electron layer is present. 
This can be realized in the 2SC, the gapless CFL (gCFL) or crystalline color superconducting phases [14,20], but not in the pure 
CFL phase. 

This paper is organized as follows. Sec. II is devoted to a qualitative discussion of photon propagation in the electrosphere. 

Here the physical characteristics of the electrosphere, the photon mean free path, the photon dispersion relation, and the role of 
the Landau-Pomcranchuk-Migdal (LPM) effect that suppresses the emission of low-energy photons are addressed. The calculation 
of the photon emissivity from the bremsstrahlung process in degenerate matter (with and without modifications from the LPM 
effect) in the electrosphere, and a characterization of the non-thermal nature of the emitted photon spectrum are detailed in 
Sec. III. This section also contains numerical results and discussion of the bremsstrahlung radiation flux, including a comparison 
with the flux of the e+e~ pair annihilation process and quark- related processes. Sec. IV contains our conclusions and outlook. 
Appendix A outlines the evaluation of an integral required in Sec. III. 



II. PHOTON PROPAGATION IN THE ELECTROSPHERE 



At the surface of the bare quark star, the density of quarks drops abruptly to zero within a layer of thickness of about 1 fm, this 
length scale being set by the range of the strong interaction. The density proflle of the electrons can be determined by solving the 
Poisson equation in the Thomas-Fermi approximation. In the plane-parallel approximation for the layer of the electrosphere, the 
electron chemical potential as a function of distance z from the quark surface is given by [21] 

Me(.) = ^^; if = A£_./|^.50i.3fl^V-. (4) 

^ 1 + V-ff /Ue(0)V2a V /"e(0) y ^' 
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where a = e^/47r is the fine structure constant. For a typical ^e(O) = (10 — 20) MeV, the electrosphere is about 10'^ fm thick or 
more. The variation of /ie across the electrosphere is shown in Fig. 1. The number density of electrons at the surface of such a 
star is rie 3 x (10^^ — lO^**) fni^'^ [17]. The electron degeneracy parameter A = nl^^h/y/nieT ^ 1 for temperatures below 10^° 
to lO^^K (~ (1 — 10) MeV), which is the case some tens of seconds after the star is born. Thereafter, electrons quickly settle into 
a degenerate Fermi sea, and their Fermi momentum far exceeds their rest mass, making it a relativistic degenerate system. 
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FIG. 1. Profiles of the electron chemical potential /ie (Eq. (4)) and the plasma frequency Up in the degenerate limit (Eq. (6)) versus 
distance z from the surface of the strange quark matter star, /xio (/i2o) is the electron chemical potential in units of 10 (20) MeV. 



The optical depth of the electron layer is determined by the mean number of scatterings of photons on electrons. Bremsstrahlung 
photons from this layer would mostly have energies much smaller than the Fermi energy (e^^ = fesTV) of the electrons. Since 

electrons can only manoeuvre within an energy window Ae ^ To^fI in a typical scattering process, and T <C Tp due to 
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strong degeneracy, the typical energy of the emitted photons would be of order Ae. The electron in the intermediate state remains 
almost on-shell. 

The ratio of the photon energy to electron mass is the small parameter that justifies an estimate of the photon mean free path 
based on the Thomson cross-section. If this ratio is of order unity or larger, the Compton cross-section should be used. For a 
given electron density, the photon mean free path is 

1 STra^ „ 

It^ , ar = = 66.54 fm^ , 5 

where ctt is the Thomson cross section. The quantity neff ~ ne{T — 0)(T/fie) is the effective number density of electrons, which 
takes into account the fact that in a degenerate gas the available scattering states are restricted to the vicinity of the Fermi surface. 
For an electron density of Ue = 3.53 x 10~^ fm~^ (corresponding to a chemical potential of /ig = 20 MeV), and T = lO^^K (8.6 
MeV), this yields a mean free path ~ 1000 fm. For higher photon energies, the mean free path using the Compton cross-section 
is somewhat larger. The electron layer, being about a 1000 fm thick, is therefore transparent to photons. 

Since the emitted photons do not scatter often enough to become thermalized, the electrosphere is a source of non-thermal low 
energy photons. We will quantify this non-thermal nature by calculating the ratio of the second to first moment of the photon 
energy spectrum and comparing it with a similar ratio (a fixed T— independent number) for the Planckian spectrum. 
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To a very good approximation, the dispersion relation of the emitted electromagnetic wave in the plasma can be taken as 
U! = + k'^Y^'^ (in c = 1 units). The plasma frequency LOp serves as a low-energy cut-off. For degenerate electrons i^T j \ie <C 1), 
the plasma frequency is determined from [22] 



2 4a 2/, 7r2r2\ 



Only photons with w > u)p can propagate without severe attenuation by the plasma. The plasma frequency, being density- 
dependent, decreases with increasing z (Fig. 1). This implies that photons emitted toward the surface from the lower layers of the 
electrosphere can escape. A large plasma frequency, however, implies an overall suppression of the emissivity. The optical depth 
of the electrosphere at any distance is determined by the plasma frequency, which in turn is a function of the electron density 
profile. The emissivity Q is thus a function of u)p{z). The total luminosity can be expressed as 

/•z=zo 

L = 4nRl Q{ivp{z))dz, (7) 

Jz=0 

where zq is the thickness of the electrosphere, and Bg the radius of the star. 

It should be noted that some (possibly paired) quark matter is admixed with electrons in the innermost region (of thickness 
about 1 fm) of the electrosphere. The screening effect of quark or electron matter will be taken into account in our calculation 
by modifying appropriately the exchanged photon that mediates the e~e~ scattering. We will see that electric screening and 
magnetic damping effects play only a small role for emission of low energy photons. Such photons can escape despite the large 
plasma frequency Up (^ 20 MeV) of the ambient quark matter because the thickness of the electron layer is about the same as 
c/iUp ^ 10^^^ cm. 

Bremsstrahlung radiation is known to be suppressed by multiple scattering of electrons (LPM effect [23]) within the formation 
time of the emitted photon. This is especially relevant for low energy photons (which have a large formation time) that are emitted 
mostly in the forward direction (small emission angles). The LPM effect is important for wave numbers k < kc that satisfy [24] 

■ (8) 



El 



PM 

where is the typical energy of the emitting particle. For degenerate electrons of energy fi^ ^ (10 — 20) MeV, the quantity 
Elpm is given by [25] 

Elpm{(N) = 3.8 X 10i2js^o(cm) . (9) 

Above, Xq is the typical radiation length in the medium. If, for simplicity, wc assume a static screened Coulomb interaction 
between electrons, the radiation length is energy independent for high energies of the emitting particles, and is given by ^ 

2 

= 4n,a37n(183) ' 

which yields £^lpm = 70 keV and kc ~ 1.4 GeV for /ig = 10 MeV. This large value of E^pm implies that the LPM effect will 
strongly suppress the number of photons radiated. For a photon of wave number k, the factor by which the bremsstrahlung process 
is reduced is given by [24] 



where in quoting the numerical result above, we have used fig = 10 MeV and k = 0.5 MeV as an example. This reduction shows 
that the inclusion of the LPM effect is important for a realistic estimate of the emissivity. The LPM suppression depends on the 
density of scatterers through /ig. 



Although we use the zero-temperature expression for Ue here, the available scattering states are properly accounted for by the Fermi 
functions in the LPM calculation to follow. 
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III. PHOTON EMISSIVITY FROM THE BREMSSTRAHLUNG PROCESS 



The photon emissivity from the bremsstrahlung process is 



^ sh 



(12) 



sptn 

X nF(wpJnF(wp2)nF(wp3)riF(wp4)(27r)^5^(Pf - Pi)S{Ef - E,) 



The subscripts i = 1 to 4 refer to electrons (1,2 label the incoming states and 3,4 the outgoing states), whereas (ti)/j,k) refer 
to the 4-niomcntum of the emitted photon. The symmetry factor ,s = 2. The phase space for electrons is convoluted with the 
appropriate Fermi distribution functions npii^pi) = l/(e'^fi/"^ + 1) (in ks = ^ units) and hp = l — np, respectively. For simplicity, 
and to highlight the role of the LPM effect, we first obtain the analytic T and /Ug dependences of the emissivity without the LPM 
suppression factor. The effects of the LPM effect are investigated in the second stage. In this case, an analytical analysis becomes 
unwieldy and hence we present results of numerical calculations. A comparison of the respective emissivities enables an overall 
assessment of the importance of in-medium multiple scattering effects. 



A. The Bremsstrahlung Process in Degenerate Matter 



For conditions of interest here, electrons in the electrosphere are strongly degenerate. The emissivity will thus be dominated 
by low energy photons (ujp < Wk <^ He)- The contribution from high energy photons (wfc me), which requires a full quantum 
calculation of the bremsstrahlung process, is expected to be small. We utilize the quasiclassical approximation for low energy 
photons, in which case the calculation of the emissivity simplifies if we use Low's theorem for bremsstrahlung [26]. 

Low's theorem states that the first two terms in the differential cross-section of bremsstrahlung, expressed as a power series in 
the energy loss Uk, is unique and exactly calculable in terms of the corresponding elastic amplitude. The presence of the factor 
ojk in the numerator of Eq. (12) ensures that in the low energy limit, the dominant contribution to Q is non-singular and can be 
estimated once the leading term in the series (~ l/uJk) is established. Retaining only the first term in Low's soft-photon expansion 
shows that the full scattering amplitude for bremsstrahlung factorizes into an elastic part, which is M0ller scattering, and a part 
which gives the classical intensity of photon emission divided by the energy Uk- At low energy, the probability of successive 
photon emission is enhanced, which causes a naive perturbative expansion in powers of the electromagnetic coupling to fail. The 
differential cross-section can alternatively be understood in terms of the mean number of photons emitted, which is equivalent to 
the classical intensity divided by the mean energy of one photon. This factorization of amplitudes simplifies our calculation to a 
large extent. The energy exchange in e~e~ scattering being small, we can write the matrix element explicitly as 



iM = Me. 



P4.k p2.k ps.k pi.k 



(13) 



where Mei denotes the interaction for elastic scattering of the incoming electrons, e'^ denotes the polarization vector of the photon, 
and e is the charge of the electron. The term in the second bracket in Eq. (13) is the familiar result from classical radiation theory 
for an electron that accelerates under an impulsive force. In the low-energy limit, the emission from the exchange diagrams simply 
doubles the overall result for the emissivity from the direct diagrams (no interference). Since there is a symmetry factor of s = 2 
in the denominator of the emissivity expression Eq. (12), we may ignore the identity of the particles altogether and evaluate the 
emissivity from the direct diagrams alone. The spin sum can be performed separately over the electrons and the photon thanks 
to the decoupling at low energy. We are thus led to the following expression for the emissivity 



^ = nj (2^(2^"^^'"-)"^^'"-^ '^''''[jJ-k 

3 E 



dflk 



A=l,2 



P4.k p2.k p^.k pi.k 



(14) 
(15) 



(dae;) = ^^(n |Y^)El^ed'(27r)3'5'(Pf-Pi)27r5(i;/-£0"F(a;p3)nF(a;,J . 

The change of variable from k to introduces the square root factor as the Jacobian, and the integration over LOk runs from 
to oo. Medium modifications to bremsstrahlung emission can be encapsulated in the refractive index of the electron plasma 

n = ^1 — w^/o;^. With the definitions p = Ev and k = w^n = Wfenn, where |n| = 1, the angular integral in Eq. (15) is 



5 



J (27r)3 



A=l,2 



E4h- 



P4 ■ ( 

1- J^v^.n 



l>2 ■ ( ' 
1 - ^ V2 



n 



p.; ■ ( 

n. "J 



Pi - (' 

1 - 42- vi 



n 



(16) 



where h± = 1 ± i and wc have neglected terms hke 7 • e* 7 • fc in the numerators compared to e* • p. Since ujk — ^ E /i,. , 
we may approximate = = 1. Although the refractive index n <C 1, we cannot neglect the factors [1 — Vj • n] in the 
denominator, since that would make the matrix element vanish by momentum conservation (^^Vi = 0). Thus, for small n, the 
leading contribution to the squared matrix element for bremsstrahlung is 0{n'^) and to the cross section O(n^), one power of n 
coming from the Jacobian in Eq. (15). This allows us to perform an expansion of the angular integral in even powers of n. To 
C(n^), and in the ultra-relativistic limit, the angular integral over dSlk can be evaluated, and one finds 



da \ _ 



= {da el) 



8a n^sin^OcM 



5n 



(17) 



where a = e^/(47r) is the fine structure constant and 6cm is the angle in the centre of mass (henceforth CM) frame. In comparison, 
the vacuum result is [27] 



d<y\ ..J F{0 

- — = {daei)\a 

dUk J 1 '^k 



(18) 



with 



m = - 



, E . (6cM 
me V 2 



The function F(^) reduces to S^^/Stt in the limit ^ <C 1 (small angle scattering in the CM frame) which is the region of interest 
here. A comparison of the vacuum and medium expressions reveals a suppression factor of order x (me/ due to the refractive 
index n of the medium. 

In general, there are two conditions that must be satisfied for the validity of the two preceding equations. The first is that 

m^l E <^ 1, because electromagnetic radiation from fast-moving particles is emitted roughly in a very narrow cone of this opening 
angle (this is equivalent to ^ ^ 1). Secondly, the exchanged photon 3- momentum q should be much larger than the change in q 
induced by soft photon emission [27]. Explicitly, 



< 0CM < I -g- 



(19) 



in the ultra-relativistic case. These limits on 0cm suffice to give results for da/duk that are correct to logarithmic accuracy. To be 

more acciirate, one needs to go beyond the quasi-classical approximation which is not required for our purposes, since we are only 
interested in the emission of low energy photons. Emission of photons with large energies is severely supressed by Pauli blocking 
in the degenerate electron sea, and makes a negligible contribution to the emissivity. 

The soft photon limit also simplifies the calculation of the in- medium elastic scattering cross section (dagi)- The scattering 
occurs via the exchange of a screened photon propagator. Based on the general principle of gauge invariance, we can use Weldon's 
decomposition [28] of the amplitude 



Mel = Mpi,P3)D''''{qo, q).Mp2,Pi): MpT.pT) = eu{Pih,MPi) . 
and write it in terms of the transverse and longitudinal (temporal) parts of the one-loop self-energy as 



Mel = J^.D>'%q)J^ = 



JoJa 



JtJt 



(20) 



(21) 



where the exchange 4-momentum 9 = (go q) = Pi ~P3 = P4 ~P2- Heiselberg and Pethick [29] have shown that in the limit of small 
90) |q| «C /ie, the magnetic part of the electron current (spin flips) can be neglected, as is evident from Gordon's identity [30]. 
In this case, the electron current reduces to J/n(p™,p?"*) — > 2ep-". The squared matrix element then becomes 



\MeiY = l& e^EiEi 



1 



+ 



(1 — a;-^)cos <f) 



(22) 
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where x = (jo/|q| and (f) is the angle between the components of the velocities of the incoming particles transverse to the momentum 
transfer q. While qq is small, the momentum transfer q can be larger since the maximum angle by which electrons can scatter on 
the Fermi surface is set by the upper bound in Eq. (19). 

The longitudinal and transverse self-energies 11^ and IIt are in general complex functions of go and |q|. These self-energies can 
be obtained in the hard dense loop approximation ^. Since the exchanged photon also has small energy, we will use approximate 
forms [31] in the nearly static limit go/|q| ^ 1: 



Hi 



I'D 



and IIt = i—-r-r 
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For degenerate unpaired quark matter 



2 2 2\ 



(23) 



(24) 



where Zi = +2/3 or —1/3 is the charge of the unpaired quark species i. 

In the 2SC phase, up and down quarks of one color are unpaired. Strange quarks can be unpaired or form a spin one condensate 
which breaks rotational invariance. For the paired quarks in the 2SC phase, photon-gluon mixing can occur. The residual U{1) 
symmetry is a combination of the electric charge Q and the SU{3) hypercharge Tg. The Debye mass associated with the photon 
is given by [32] 



rrijj = 



37r2 



where Hg is the quark chemical potential, of order hundreds of MeV. For degenerate electrons, 

2 2 2 /o 



(25) 



(26) 



In what follows, we denote the phase dependent screening mass by mu, noting that its magnitude is generally larger than 
rrie- Contributions to the emissivity from screened electric effects will turn out to be much smaller than those from unscreened 
magnetostatic modes. 

It is easier to evaluate duei, a Lorcntz invariant quantity, in the CM frame, and later transform to laboratory co-ordinates 
for the final phase space integrations. In the CM frame, elastic scattering between identical particles implies that go = and 
cos (f) = —1. If the scattering angle 9iab is small, so is the scattering angle 6cm in the CM frame since 6cm = ^6iab for identical 
particles. Moreover, g^ = — = —EiEsd'^,^^. The elastic differential cross-section then simplifies to 



dae 



1 1 r d 

2E^2E^J (2^ 



d^P3 d^pi 



IQe^ElEl 



+ ■ 



E1E36CM + ml E1E36CM 



(27) 



X {2Tryd%Pi + Pa - P3 - P4)27r5(^i + E2 - E3 - Ei)nF{E3)nF{EA) 



Denoting the total CM energy as E1+E2 = 2£ and using the momentum conserving delta- function to perform the d^pi integration, 
we find 



El 



1 



1 



ErE36cM +ml, E1E36C 



M 



5{£ - E3)(hF{E3)f 



(28) 



where a = /An is the fine structure constant. Since the main contribution comes from the vicinity of the Fermi surface, we can 
change variables as d^ps « E^dEsdQp^ (the mass of the electron can be safely ignored here). The delta function facilitates the 
integration over rfjpal with the result 



da 



'P3 



1 



'CM 



+ 



1 



'cm 



a = mD/£ 



(29) 



Without loss of generality, we can take Pi and P2 to be along z and — z, respectively, so that rfOpg = dQ.cM- It can be checked 
at this stage that wc recover the vacuum result for M0ller scattering. If the blocking factors from the above expression are set to 
unity, and a = 0, we find 



^ We are ignoring the smaller contribution from finite temperature eflfects (hard thermal loops) . 
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da^i 4a^ ^^^^ 

which is the most singular part of the M0ller cross-section in the ultra-relativistic hmit. (In making the comparison, a symmetrizing 
factor of 1/2 has to be reinserted, because particles in the initial state are indistinguishable). Returning to the in- medium expression 
Eq. (29), wc note that the integral over ddcM should also include the ^cM-dependence from F{^). Integrating (j) from to 2n, 
and 0CM from Omin to 0max as prescribed by Eq. (19), wc obtain 



2 V 6'i.„ 2 \ 91 



92 . 

rmn 



(31) 



The prime on da^i is to remind us that the angic-dcpcndent part from photon emission must also be included in performing the 
angular integral in the elastic cross-section. Next, the du)k integral in Eq. (14) is to be performed. The lower limit on Wk is 

^min = while the upper limit is u>max ~ ^a;| -h ml. The upper limit is set by noting that for small angles, we can write 

l^lmai; ~ Ed„iax ~ E{m^/E) = m^. Although photons of energy LOk > L^max arc also emitted, their contribution, which cannot be 
included consistently within our approximations, is expected to be small in the degenerate limit. 
The integral over Uk can be rewritten as 



where 



/i=ln(-^) , (33) 



rrieUJk 



In '^^ ' -In 



(34) 



The Debye mass mD ^ rUe, which allows us to drop contributions from I2 and I3, and retain the largest contribution from /i. 
This contribution is easy to evaluate, since uip is a slowly varying function of the distance z for temperatures less than 1 MeV [21]. 

The logarithmic factor in Eq. (33) is only weakly dependent on Uk as compared to the square root prcfactor and can therefore be 
pulled out of the integral with only a slight underestimation of the integral. This leads to the expression 



0J„ 



Inserting this result in Eq. (14), the emissivity becomes 

64a3^e f d^pi d^p2 , , , sfiFiS)'^. f E\( 1 u^l 3wp 
Q= / 777-T?777-T?'^J'Ki)"J'K2) ci — i ^ + o . .9 , „ . ' TTZ^tan [ — ]} • (37) 



bfi J (27r)3(27r)3 £^ \me ) \ 2 u^^ + ml 2me 

We note that electric effects are not very important for the conditions of interest here, since they are strongly supressed because 
(me/mr>)^ ^ 1- We have therefore dropped contributions of Oi ml/m'jy 1 . As we have focused on soft exchanges, the magnetic 



damping for such modes vanishes; consequently they provide the largest contribution to the integral. 

For strongly degenerate electrons, further reduction of Eq. (37) is afforded by the nearly step function character of the Fermi 
functions. We can approximate £" w /ie, so that the logarithmic factor is now a constant and can be pulled out of the integral. 
The variable £ can be converted to laboratory coordinates using the relation 



COS 9i2 = COS Oicoa 82 + sin 0isiu 6*2008 {(f>i — 0; 



'2 



We can further set |Pi|, IP2I ~ Pf^, the electron Fermi momentum so that hF{£) depends only on the angular variables 9i and 
62- The integral over the momenta can be evaluated approximately in the degenerate limit: 
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(39) 



Terms involving the electron mass can be dropped because they appear only beyond second order in the above expansion. More 
details on the evaluation of the angular integral are given in Appendix A. Here, we quote the final result for the emissivity in two 
relevant limits: 



Q 



1{T, Me 



64Q:-V/i,J'(^;,)(,V(r,/(,))^ 



1 

1 + 



5n(27r)6 

, ,2 



ln(^)X(r,/.e) 



(40) 



3a;. 



2 + 



2mp 



P tan-i 



4T 



In 2 



2 ) ' 2/i, 



< T< Me 



' 2^e 



> . 



Note that ^(ojp) — » 1 for Wp <§; rrig and T{ijjp) (mc/ujp)'^ for ojp ^ rrie- 

Equation (40) provides the photon emissivity to logarithmic accuracy. A better estimate can be obtained by a numerical 
integration of Eq. (37) which can be reduced to a five-dimensional integral by choosing, for example, the azimuthal angle 02 = 0. 
The limits on then run from to 27r, and 6\ and 62 each range from to tt. In the degenerate limit, the Fermi functions 
cause the integrand in Eq. (37) to be strongly peaked for radial momenta pi ~ J52 ~ ^^e■ In this case, it is advantageous to use an 
adaptive Monte-carlo method such as VEGAS [33], which is based on stratified importance sampling. We can restrict the range 
of integration for the radial momenta from /Xg — (5 to /Ug + ^, where 5 ^ He'is chosen to be a suitable multiple of the temperature. 
This procedure provides an accuracy of better than a few percent. 

If the temperature becomes comparable to /Kg, electrons cannot be considered as degenerate and our approximations would not 
be valid. In our numerical calculations therefore, temperatures are chosen to be much smaller than /Xe- Within the temperature 
limits that set the validity of our expressions, one can compare the photon emissivity from e~e~ bremsstrahlung to that of e+e" 
pairs [34]. Such comparisons are presented below. 



B. Comparison with Other Processes 



For purposes of comparison and easy estimation, it is useful to express all fluxes in terms of typical temperatures and chemical 
potentials. We will employ Tg to denote the temperature in units of lO^K, and fi\o to represent the electron chemical potential in 
units of 10 MeV. The result for the bremsstrahlung radiation flux from Eq. (40) can then be recast as 



Jz=0 



dzQ; g = 1.14 X 10'*2<1 1-t- 



X ln(19.6Mio) X 



0.15m 



2 
10 



0.3mio + 0-26 

[ ln2 - 0.5 



— l.Gl/xiotan 



0.02 e-o-!^/'^^''"'; Tq/mio « 0.15 




(41) 



where zq is the thickness of the electrosphere. The plasma frequency of a highly degenerate electron plasma depends linearly on 
the chemical potential (see Eq. (6)) and has been rcscalcd accordingly (LUp{z) = 0.55fiio{z) MeV). Excepting the outermost regions 
of the electrosphere, we can safely assume the electrons to be degenerate for all chemical potentials and temperatures considered 
in this work. The flux (similarly re-expressed) from the pair creation process is [34] 



F± 



rZ=Zo 

/ dz Q± ; 

Jz=0 



Q± = 4.52 X 10*^{0.51 + 0.0862T9}mio79 

X exp 



1 + 



256 Vmio 



(42) 



11.83\ ^/ 11.6mio ^ _3 _i 



where 



^, , lx3ln(l + 2/a; 



3 (1 -h 0.074a;)3 6 (13.9 + xf 



(43) 
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The blackbody fliix is given by the Stefan-Boltzmann law: 

Fbb = 5.67 X lO^^Tg* erg cm-^s"^ . (44) 

For illustration, wo choose Tg = 1, an electrosphere thickness zq = 1000 fin and /iio(O) = 1. With the electron density profile 
given by Eq. (4), we find 

F = 9.06 X 10^^ erg cm-^ s"^ (45) 
F± = 1.45 X 10^® erg cm-^ s"^ 
Fbb = 5.67 x 10^^ erg cm^^ g-i ^ 

These numbers demonstrate that the bremsstrahlung radiation flux from e~e~ scattering in the electron layer dominates over that 
from pair creation. 

In Fig. 2, wc compare the flux from the brcmsstrahhmg process to that from the pair creation process, both scaled to the 
blackbody limit. For simplicity and for the sake of comparison with the analytical result Eq. (40), we assume a uniform chemical 
potential. In this case, the integral over z gives simply Zq, the thickness of the electrosphere. To be consistent, we take the 
thickness to be zq = 100 fm for which the chemical potential changes only by about 20% . The analytical estimate in Eq. (40) 
agrees well with results from numerical integrations of Eq. (37). For Tg < 2, the flux from the bremsstrahlung process is more 
than that from the pair creation process (it may even dominate at somewhat higher temperatures once the contribution from 
higher energy photons is included.) whereas for temperatures Tg > 2, the flux from the pair creation process exceeds that from 
bremsstrahlung. (This conclusion is not significantly modified upon the inclusion of LPM effects, as we will see.) Thus, for electron 
plasmas below a temperature of ^ 0.1 McV, brcmsstrahhmg radiation is an important channel for photon cooling of bare quark 
stars. Its dependence on the third power of the electron chemical potential (or linear in electron density) distinguishes it from 
other processes. 

With increasing electron degeneracy (value of He/T) or electrosphere thickness zq, the flux from brcmsstrahhmg rapidly increases 
until it starts to become comparable to the blackbody limit. At this stage, the inverse process, photon absorption on electrons 
and rescattering (Comptonization) becomes important and the opacity of the layer increases. Detailed balance then implies that 
the blackbody limit is saturated, and the electrosphere radiates like a blackbody. The inverse processes have not been taken into 
account in the curves of Fig. 2, hence they appear to exceed the blackbody limit for some temperatures. However, there exists a 
window, 0.1 < Tg < 1, in which the inverse processes can be neglected and in which the bremsstrahlung flux exceeds that from 
pair creation. 



C. Spectral Features and Typical Energies 
The spectrum of low energy bremsstrahlung photons can be read off from Eq. (32). Explicitly, 



where in the term involving the logarithm the electron energy E has been replaced by fie- In the quasiclassical approximation 
employed here, support for this spectrum exists only between uJmin = and u>max = {^p + tnlY^^- One must go beyond the 
quasiclassical approximation in order to obtain the spectrum for w > oJmax- In vacuum, the low energy bremsstrahlung spectrum 
for frequencies — tJk ^ rUg reads [27] 



1/2 

x = l-^. (47) 

Me 



For typical values of Me encountered across the electrosphere. Fig. 3 shows the role of medium effects by comparing the low energy 
spectrum with the low energy spectrum neglecting all medium effects (except the plasma cutoff). 

Although photons are emitted from thermal electrons, the non-thermal character of the bremsstrahlung radiation is evident from 
this figure. While this is a natural feature of any spectrum calculated with a low-energy cut-off, the large mean free path of photons 
in the electrosphere, and medium modifications reinforce the same physical conclusion even if the entire range of frequencies is 
covered. The extra factor of induced by bremsstrahlung in the medium produces a harder spectrum than in vacuum, while 
the total flux is smaller. The spectrum of high energy photons is expected to fall off steeply due to the large electron degeneracy. 
However, as mentioned previously, this feature cannot be calculated within the quasiclassical approximation made here. 

We turn now to estimate the characteristic energies of the emitted photons, and to contrast them with those characteristic of 
thermal blackbody radiation. In order to perform these tasks, we examine the quantities 
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FIG. 2. Photon fluxes from the bremsstrahlung and pair creation processes scaled to the blackbody limit as a function of temperature in 
units of lO^K (Tg). The curve marked (num.) shows results of numerical integrations of Eq. (37) which is to be compared with the results 
of the analytical expressions (anl.) from Eq. (40) for low and high temperatures. The (uniform) electron chemical potential is ~ 10 MeV 
with zo = 100 fm. 




(2^(2^3' 



(48) 



where F is the number of bremsstrahlung reactions per unit volume per unit time (the rate) , Q is the energy loss in photons per 
unit volume per unit time (the emissivity), and is the squared energy loss per unit volume per unit time. Good measures of 
the average photon energy and the average squared photon energy are provided by the ratios 



{u) = Q/F and ( 
It is also instructive to examine the dimensionless ratio 



= QVr. 



(49) 



(50) 



which can be used to gauge the gross spectral features of the energy spectrum. The quantity R takes into account the full effects 
of the transition amplitudes and 4-momentum conservation. 

The calculations of F (this is finite because of the finite plasma frequency) and for the bremsstrahlung process can be 
performed following the procedure adopted for the calculation of Q in Sec. Ill A. The results are 



H = 



Q = 



±ln 



t2 2nie 



+ Ami 

3(l^2+ ^2)3/2 



(51) 
(52) 



where the quantities Q and J-' are given in Eq. (40). The average photon energy is therefore 
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FIG. 3. The energy spectrum of bremsstrahlung photons at the indicated values of the electron chemical potential /Xe- The full curves are 
the spectrum from Eq. (32), while the truncated curves correspond to Eq. (47) for the low energy spectrum in vacuum valid for frequencies 
satisfying the condition — cok me- 

(lu) = — ^ I ^ "^^ (53) 

^ ' 7Y \me/\n{2me/LLJp) for < me ' ^ ' 

In Fig. 4, we show the mean energy of photons versus distance in the electrosphere for /ie(0) — 10 MeV and T = 1 MeV. The solid 
curve shows the result obtained by using the analytical expression for (a;) = J^/Ti.. Results obtained by numerical integrations 
of Eq. (37) and its counterpart for T are shown by the dashed curve. It is clear that the analytical results give a fairly accurate 
representation of the numerical results for values of Wp sampled across the electrosphere (see Fig. 1). The mean energy of photons 
varies from 0.4 MeV to 0.7 MeV, depending on the depth from which the photons are emitted. On average, this is larger than 
the mean energy of 0.5 MeV from the annihilation of e+e^ pairs. This is also very different from the spectrum of photons from a 
normal cooling neutron star with 0.1 < E/keV < 2.5. 
The ratio R is given by 

gji j 1 for LOp > me 

'^=^-^\ihr(^) forc.p«me (^4) 

In order to highlight the differences between the bremsstrahlung and blackbody radiation, the quantities {to) and R can be 
contrasted with the ideal gas benchmark quantities. 

«,.i;^lM ) Jo I" ^ ,55) 







where is the equilibrium Bose-Einstcin distribution function characterizing in-medium photons with mass ujp and zero chemical 
potential. Explicitly, 

Up + |T for Wp > T 

r(4)C(4) ^ „ _ ^ , ^ (56) 
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u(0)=10MeV, T=l MeV 
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FIG. 4. The mean photon energy as a function of distance in the electrosphere. The sohd curve shows results from the analytical expression 
in Eq. (52), while the dashed curve is obtained by numerical integrations of Q in Eq. (37) and its counterpart P. Numerical values of the 
electron chemical potential fj,e and temperature T are as shown in the inset. 



+ 3rnT 



for LOr,^ T 



r(5)C(5) 
r(3)c(3) 



~ 10.35 r2 for w„ < T 



These results imply that [35] 



Rid 



1 for Ldp->T 
1.42 for cj„ < T 



(57) 



(58) 



The behaviors of R versus tOp/nif, and Rid versus LOp/T are compared in Fig. 5, which highlights the non-thermal nature of the 
bremsstrahlung radiation emitted from degenerate electrons. In light of the limited frequency support for bremsstrahlung, a more 
rigorous quantum calculation is called for to characterize the spectrum more accurately. 

D. Inclusion of the Landau-Pomeranchuk-Migdal Effect 



As noted in Sec. II, multiple scattering of electrons within the formation time of the radiated photon can lead to the so- 
called LPM suppression of bremsstrahlung. We had omitted a detailed consideration of this effect in order to isolate the T and 
/ie-dependences of the emissivity in the absence of the LPM effect. Since an analytical estimate including LPM suppression is 
cumbersome, we have performed a numerical integration of Eq. (37) utilizing the approximate suppression factor SLPuik) from 
Eq. (12). 

In Fig. 6, the effects of the LPM effect are assessed. The LPM effect reduces the simple analytical estimate in Eq. (40) by 
a significant amount. Specifically, for /iio — 1, the reduction factor (relative to the emissivity without the LPM effect) is ~ 60, 
while for /iio = 2, the reduction factor is ^ 350. The temperature range is 0.1 < Tg < 10. The suppression is stronger at higher 
densities (compare the left and right panels), since the density of scatterers controls the LPM effect. The suppression factor for 
/iio = 0.5 (not shown in the figure) is approximately a factor of 10. 

The mean photon energy Q/F does not change much from the analytical estimate upon the inclusion of the LPM suppression 
factor. This is because both the rate and the energy loss are affected similarly and therefore the LPM effect nearly cancels in the 
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FIG. 5. A comparison of the ratio R in Eq. (54) for bremsstrahlung photons and Rid in Eq. (55) for thermal photons with in-medium 
photon mass uip. 

ratio. For the same reason, the LPM suppression does not significantly change the shape of the spectra discussed in the previous 
section. The spectra are controlled primarily by the ratio of plasma frequency to electron mass. 

Equilibrium radiation from quark matter and bremsstrahlung radiation from quark-quark collisions in the surface layers also 
produce photons. In Ref. [36], luminosities from these processes were compared to the blackbody luminosity; it was found that 
for T < 2 X 10^*^ K (i.e., Tg < 20), the equilibrium radiation is negligible. For these temperatures, the non-equilibrium radiation 
from quark-quark Bremstrahlung dominates, although it is still 4 orders of magnitude less than the blackbody luminosity (even 
this may be an overestimate according to Ref. [34]). Accounting for the multiple scattering (LPM effect) of quarks and absorption 
by the electron layer leads to the conclusion that an additional suppression by 2 orders of magnitude is likely [37] , so that overall, 
the photon emissivity from quark matter is about 6 orders of magnitude smaller than the blackbody emissivity. The dominant 
contribution thus comes from e"*"e~ pair production via the Schwinger mechanism in the presence of strong electric fields or the 
e~e~ bremsstrahlung process in the thin electron layer. The relative importance of these two processes depends sensitively on the 
temperature. 

In the range T < lO^K, bremsstrahlung is the dominant emission process, since the emissivity from the pair creation process is 
vanishingly small. The analytical treatment of the bremsstrahlung process performed in this work provides a good approximation 
to the numerical results in this range of temperatures. At low temperatures (T < 10^ K), the bremsstrahlung rate is exponentially 
small as seen from Eq. (40). For temperatures in the range T ^ 10*K, the inverse process of photon absorption also becomes 
important as the rate approaches the blackbody limit. The main conclusion is that for temperatures T < lO^K, the bremsstrahlung 
process is likely to be the dominant photon emission source, even in the presence of LPM effects, when compared to the annihilation 
of e+e^ pairs. 

IV. CONCLUSIONS AND OUTLOOK 

We have calculated the photon emissivity from the e~e~ e~e~7 bremsstrahlung process occurring in the electrosphere of 
a bare strange quark star. For degenerate electrons and for low energy photons, unscreened magnetic interactions provide the 
largest contribution to the emissivity, while the contribution from screened electric interactions is comparatively negligible. For 
temperatures in the range 10®K< T < lO^K, the emissivity from the bremsstrahlung process surpasses that of the e+e" pair 



14 




FIG. 6. Photon fluxes from the bremsstrahlung and pair creation processes scaled to the blackbody value as a function of temperature 
in units of lO^K (Tg). In the case of the bremsstrahlung process, the upper (lower) curve shows results without (with) the effects of the 
Landau-Pomeranchuk-Migdal effect. Numerical values of the electron chemical potentials are as indicated. 



production process occurring in the presence of a strong electric field in the electrosphere. Within the soft photon approximation, 
the mean energy of the bremsstrahlung photons escaping the electrosphere is ~ 0.5 MeV. This estimate may be slightly altered if 
the entire spectrum is sampled. Under strongly degenerate conditions, the contribution of high energy photons to the emissivity 
is expected to be small in comparison to that from low energy photons. 

Photons from the bremsstrahlung process are non-thermal for sufficiently thin electrospheres and moderate electron chemical 
potentials, because such photons do not undergo sufficient number of scatterings in the electrosphere. Inclusion of the Landau- 
Pomeranchuk-Migdal effect, which arises due to multiple scatterings of the emitting electron, reveals that even in the presence of a 
large suppression from this effect the bremsstrahlung process emerges as the dominant source of radiation from the star's surface 
for temperatures in the range 10®K< T < lO^K. We may regard the LPM-inclusive bremstrahlung rates as a conservative estimate, 
since we have approximated photon exchange by a static Coulomb field. If retardation effects were included, the suppression will be 
smaller leading to somwewhat higher rates. Both the "bare" and LPM-corrected rates take into account the dielectric suppression 
effect of the electrosphere. The consequence is that the LPM suppression is two orders of magnitude less than that in the case the 
photon is emitted into vacuum. There is an additional mitigating effect of about 15% on the LPM suppression due to the strong 
external electric field, which can bend the electron trajectory causing the formation time of the emitted photon wave packet to 
decrease. This implies that the emission angles can be increased beyond 9max = m^/ E. However, this only changes the emissivity 
by a pre-factor ~ 4/3, and to logarithmic accuracy our principal results are valid even in the presence of strong electric fields. In 
view of this, our LPM-corrected rates should be treated as a lower bound. 

An intriguing result of our calculations is that for T < 10~^ /^e(0), bremsstrahlung emission becomes comparable to blackbody 
emission. This indicates that bremstrahlung absorption has to be taken into account in estimating the optical thickness of the 
electrosphere (see the discussion following Eq. (5)). This will lead to better estimates of the emissivity applicable in the temperature 
window excluded by Eq. (42). The remaining task is to perform calculations of the thermal evolution of strange quark matter 
stars including the non-thermal bremsstrahlung and thermal pair annihilation processes in order to provide baseline calculations 
of the surface luminosity versus age. Results of such calculations will be reported separately. 
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APPENDIX A: ANGULAR INTEGRATION 



The factor X(T, /Xg) in Eq. (40) is evaluated below. Since only the relative angle between 6i and 62 enters in the integration, and 
we integrate over all relative orientations, we can choose a direction pi = z to decouple the two integrations. Then, 



dOi / (i(cos 



£2 



(Al) 



dx 



Jj^Q{He-l^{l-x)+ml]^/^)/T^ 



-{1 — x) + ml 



Since Pf^ ^ and T, the dominant contribution comes from the region aromid x « — 1; elsewhere the integrand is exponentially 
suppressed. Therefore, we can restrict the range of x from -1 to —1 + 6, where 6 is small. Defining e = (1 +x)/2, we can integrate 
e from up to a characteristic 5 ^ 10(r//Xe) when the exponential suppression becomes effective. This leads to 



167r^ 



de 



1 



(1-e) \l + et'''^/^T 



ln2 - - 
2 



(A2) 



This result is valid as long as mg/2neT < 1 and is accurate up to corrections of higher powers of T/^e- Therefore, it is a 
better approximation for large electron degeneracies. When m\/2iieT 1, the blocking factors approach zero due to increasing 
degeneracy of the Fermi sphere, and we obtain 



(A3) 



This result shows that the emissivity is exponentially small at very low temperatures. 
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